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ABSTRACT 

We predict the evolution of giant clumps undergoing star-driven outflows in high-z 
gravitationally unstable disc galaxies. We find that the mass loss is expected to occur 
through a steady wind over many tens of free-fall times its ~ 10 Myr) rather than 
by an explosive disruption in one or a few tg. Our analysis is based on the finding 
from simulations that radiation trapp ing is negligible because it destabilizes the wind 
(jKrumholz fc ThompsoJ l2012l 120131 ). Each photon can therefore contribute to the 
wind momentum only once, so the radiative force is limited to L/c. When combining 
radiation, protostellar and main-sequence winds, and supernovae, we estimate the 
total direct injection rate of momentum into the outflow to be 2.5 L/c. The adiabatic 
phase of supernovae and main-sequence winds can double this rate. The resulting 
outflow mass-loading factor is of order unity, and if the clumps were to deplete their 
gas the timescale would have been a few disc orbital times, to end with half the 
original clump mass in stars. However, the clump migration time to the disc centre 
is on the order of an orbital time, about 250 Myr, so the clumps are expected to 
complete their migration prior to depletion. Furthermore, the clumps are expected to 
double their mass in a disc orbital time by accretion from the disc and clump-clump 
mergers, so their mass actually grows in time and with decreasing radius. From the 
6-7 giant clumps with observed outflows, 5 are consistent with these predictions, and 
one has a much higher mass-loading factor and momentum injection rate. The latter 
either indicates that the estimated outflow is an overestimate (within the \-a error), 
that the SFR has dropped since the time when the outflow was launched, or that the 
driving mechanism is different, e.g. supernova feedback in a cavity generated by the 
other feedbacks. 
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1 INTRODUCTION 

In our developing picture of violent disc instability 
(VDI) at high redshift, the gas-rich discs fed by 
cosmological streams give birth to giant baryonic 
clumps that are the sites for intense star forma- 
tion. The clumps are expected to migrate toward 
the disc centre on an orbital tim e scale where 
they coalesce into the central bu lge (iNoguchil 19991: 
Immeli et al. 20041: Bournaud. Elm cgrccn fc ElmegreenI 
2007t lElmegreeirBournaud fc Elnicgree^ 120081: 
Genzel et all 120081: iDekel. Sari fc Ceverind 



EQOI; 

Ceverino. Dekel fc Bournaud ^1(|). This was pro 



20091: 



Agertz. Tevss icr fc Moor^ 



posed as a mechanism for the formation of galactic 
spheroids, in parallel with the traditi onal scenario 
of s ph eroid formation by merger s (jGenzel et al.l 
120081: IDekel. Sari fc Ceverind l2009t) as well as 
a scenario for the formation of g lobula r clusters 
([Shapiro. Genzel fc Forster Schreiber 201(il). and for 
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feedin g the central black holes (jBournaud et al.|[20Tll 
l2012f ). However, stellar feedback can generate outflows 
from the clum ps (Mur ray. Quataert fc Thompsonll2010l: 
IKrumholz fcDekel,2010f K These outflows were assumed 
to be very intense on a timescale of a few free-fall times 
and thus lead to signi ficant mass loss and possibly to 
clump disruption (jMurrav. Quataert fc ThompsonI 
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20101: iHopkins et all 120121: iGenel et all I2012D . 
Murray. Quataert fc Thompsonl ( 2010[ ) argued that 



the high-z giant clumps are hkely to be disrupted 
by momentum-driven feedback, as are their smaller 
counterpart molec ular clouds in the Milky Way at 
low redshift, but iKrumholz fc Dekell (|2010f ) pointed 
out that this would be possible only if the efficiency 
of star-formation per free-fall time es is significantly 
higher than the value implied by observations of both 
nearby and high-z galaxies (IKrumholz fc TanI 120071 : 
IKrumholz. Dekel fc McKei |2012|) . namelv the value 
associated with the Kennicutt- Schmidt (KS) relation. 

iGenzel et all (l201lD reported pioneering observa- 
tional evidence for outflows from giant clumps in five 
z ~ 2 galaxies. The SFR is estimated from the 
He luminosity. The clump properties of radius and 
characteristic velocity are measured directly, or al- 
ternatively the gas mass is derived from the SFR 
assuming that the KS rela tion observed on sub- 
galactic scales at z = P (IKrumholz fc TanI l2007l : 



IKrumholz. Dekel fc McKeil2012 ) and on galact i c scales 
at z ^ 2 (IGenzel et all 120101: iDaddi et all l2010bl: 



Tacconi et al.H2010l : iDaddi et al.l FoiOat iTacconi et al l 
2013[ ) also holds within ind ividual giant clurnps at 



(see preliminary results bv iFreundlich et 311120131 ). The 
outflow velocity is evaluated based on both the centroid 
blue-shift and the width of the broad-line component 
of the Hq, emission, and the mass outflow rate is esti- 
mated using a vari ety of alterna t ive m odels. Based on 
these observations, IGenzel et al.l (j2011[) estimated that 
the typical clumps in their sample drive winds of mass 
loading factor ~ 1, namely a mass outflow rate that 
is comparable to the SFR. These winds arc expected 
to deplete the clump gas in several hundred Myr af- 
ter turning about half the clump mass into stars. Two 
extreme cases, both in the same galaxy, ZC406690, in- 
dicate stronger outflows, with 77 ^ 3 — 7 and depletion 
times of (100 - 200) Myr with only 12-25% of the origi- 
nal clump mass turning into stars. In the typical clumps, 
the actual driving force of the outflows is (2 — 4)L/c, 
where L/c is, the contribution of a single-scatter radia- 
tion force, while in the most extreme outflow it is es- 
timated to be as large as 34 L/c (though with a very 
large uncertainty). 

iNewman et al.l ([2012!) performed follow-up obser- 
vations at higher resolution on the galaxy with the ex- 
treme clumps, ZC406690. They compared two clumps 
in this galaxy, both of which are driving winds, but 
with very different properties. One of the clumps shows 
a considerably larger mass, energy, and momentum flux 
than the other. They propose that these two clumps 
represent different evolutionary stages of the same phe- 
nomenon, and that the more energetic of the two out- 
flows cannot easily be explained by any of the wind 
launching mechanisms that have been proposed in the 
literature. 

In this paper we seek to provide a unified frame- 
work for comparing different potential outflow launch- 
ing mechanisms, and then use this framework to predict 
the outflow properties expected from stellar feedback, 
and understand what can be learned from the obser- 



vations conducted to date. We consider the momen- 
tum injected into the wind by momentum-conserving 
stellar feedback mechanisms, and by the more energy- 
conserving supernova feedback. We predict the expected 
mass loading factor and momentum injection efficiency. 
We focus in particular on the question of whether the 
migrating clumps arrive at the centre massive and in- 
tact or lose most of their mass to outflows while still in 
the disc. 

In we develop a simple theoretical framework 
for dealing with the momentum that drives outflows 
from star-forming clumps. We present in comparison 
the timescale for clump migration and the accretion rate 
into the clumps during migration. In fJ3]we go through 
the momentum budget for the outflows. In fj4l ad- 
dress the implications for the evolution of high-z giant 
clumps. In 3SI we compare the observational estimates 
to the predictions, and discuss the implications on the 
outflow driving mechanisms. In fJSjwe conclude our re- 
sults and discuss them. 



2 THEORETICAL FRAMEWORK 

2.1 Momentum versus Energy Feedback 

Newly formed stars generate outflows by injecting mo- 
mentum and energy into the interstellar gas. Our goal in 
this section is to develop a basic theoretical machinery 
to describe this phenomenon. We flrst address the roles 
of momentum and energy in this context, and clarify the 
terminology of momentum-driven versus energy-driven 
feedback. 

The central conceptual challenge is that cool in- 
terstellar gas is highly dissipative, so energy is always 
lost to radiative processes. Indeed, in cold gas the cool- 
ing time is almost always short compared to dynamical 
timescales. Thus in launching a wind, what we are in 
the end always concerned with is the amount of mo- 
mentum that is transferred by stellar feedback to the 
gas [3 We can address two extreme ways for this trans- 
fer to occur. First, momentum- conserving (PC) trans- 
fer, where ejecta from stars (photons, winds, supernova 
ejecta) collide with the ISM inelastically, transferring 
its momentum but losing some of its energy. Second, 
energy- conserving (EC) transfer, where stars heat the 
interstellar material, either radiatively or via shocks, 
to temperatures high enough that the cooling time be- 
comes much longer than the dynamical time. When this 
happens the hot gas expands adiabatically, transfer- 
ring momentum to the cool phases of the ISM as it 
does so. We generally refer to the former mechanism 
for launching a wind as momentum-driven and the lat- 
ter as energy-driven, but this is a somewhat mislead- 
ing nomenclature, because the rapid cooling in the cold 
phases of the ISM implies that in either case what ul- 
timately matters is the momentum transferred to the 
cold gas. 



^ Note that when considering a spherical outflowing shell one 
refers to the overall momentum in the radial direction, which is 
not necessarily conserved. 
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The energy-conserving case is in general much more 
efficient. To see this, consider a source of ejecta with 
outflow velocity Vg (T4 = c for radiation) and mass flow 
rate Ms (or the equivalent energy outflow rate for ra- 
diation). In the PC case, after time t, the ejecta has 
pushed a wind of mass Mp and velocity Vp obeying 



MpVp = MstVs . 



(1) 

In the EC case, the ejecta has pushed a wind of mass 
Me and velocity Vg obeying 



(2) 

The ratio of wind energies between the EC and PC cases 



El 



MpVp^ 



MA 



(3) 



and, more importantly, the corresponding ratio of mo- 
menta is 



MeK 

MpVp 



1/2 



> 1. 



(4) 



Both ratios arc much larger than unity as long as the 
wind mass is much larger than the mass of the direct 
ejecta from the source, and the wind velocity is much 
smaller than the original ejecta velocity. From eq. ^ 
we learn that the efhciency of injecting momentum into 
the wind, ip = MV/ A/IstVs, is much larger in the EC 
case than in the PC case. It turns out that radiative- 
pressure and stellar winds are PC mechanisms, while 
the expanding supernova ejecta have an early adiabatic 
phase that makes them closer to EC. 

It is worth pausing for a moment to elaborate on 
how it is possible for the radial momentum imparted to 
the wind to greatly exceed that provided by the source 
when the flow is energy-conserving, and what distin- 
guishes the EC and PC cases. The characteristic sig- 
nature of the EC case is the presence of some mecha- 
nism that carries information and forces between differ- 
ent parts of the expanding shell of swept-up gas, allow- 
ing them to push off one another and thereby greatly 
increase their radial momenta while leaving the vector 
momentum of the shell as zero. The mechanism respon- 
sible may be sound waves traveling through hot gas that 
communicate forces via gas pressure, it may be pho- 
tons bouncing from one side of a spherical shell to the 
other that carry information via radiation pressure, or it 
may be something else, as long as that something allows 
forces to be transmitted from one side of the expand- 
ing shell to the other. In contrast, the distinguishing 
feature of PC flows is that there is no causal commu- 
nication between different parts of the expanding shell, 
and as a result they cannot push off each other and 
increase their radial momenta. 



2.2 Momentum Injected by Star Formation 

Star formation in clumps, either giant clumps at z ^ 2 
or much smaller star-forming clumps in the nearby uni- 
verse, generally does not continue until it has exhausted 



the available gas supply into stars. Instead, it ends after 
some fraction of the initial gas has been converted to 
stars. 



M, 
M, 



(5) 



The remaining gas is either directly expelled by the 
winds, or is removed by tidal stripping or simply drifts 
off after the winds have eroded the clump mass enough 
to render it unbound. The quantity e^, is commonly re- 
ferred to as the star formation efficiency. The time of 
final gas depletion is tdep = M^/Mi,, where M^ is the 
star formation rate (SFR) and M^ is the final stellar 
mass at tdep- This expulsion can occur via a gradual 
wind that removes mass from the clump continuously 
as it forms stars, via an explosive event that removes 
the bulk of the mass on a timescale comparable to or 
smaller than the clump dynamical time, or some combi- 
nation of the two. Our goal in this section is to develop 
some basic theoretical machinery to describe this phe- 
nomenon. 

Let Mc and Rc be the mass and radius of a star- 
forming clump, and let 



iff 



Rc 



GMc 



GAL 



Rc 



(6) 



be the corresponding free-fall time and characteristic 
velocity, respectively, where we are dropping factors of 
order unity for simplicity^ The crossing time Rc/Vc is 
of the same order as ts, and the escape speed is of the 
same order as 140 The instantaneous SFR is written as 



Mc 
iff 



(7) 



(8) 



The corresponding SFR timescale is 

_ Mc _i 

tsh- = —7— — Eff tg . 
M* 

The quantity eg is known as the SFR efficiency per 
free-fall time, or the rate efficiency, to distinguish it 
from the overall star formation efficiency e*. Obser- 
vations of star formation in a wide variety of en- 
vironments at a wide variet y of redshifts strongly 
constrain that eg ^ 0.01 (Krumholz fc TanI 120071: 
iKrumholz. Dekel fc McKeell2012[ ). In particular, there 
are indications from observed CO that th e same rela- 
tion with the same eg is also valid at z ~ 2 (jPaddi et all 



^ The expression tg = ^ "iix j [?>2G p) , where p is the mean density 
within the clump, is \/S/-k times the expression in eq. ((gJ. 
^ For an object with zero pressure and magnetic field the escape 
speed is \/2Vi, but observed molecular clouds in the local uni- 
verse are only magnetical ly supercritical by factors of ~ 2 (e.g. 
iTroland fc CrutcheH [20o3 ) . corresponding to a reduction in the 
escape speed to Vc- There are no direct measurements of mag- 
netic field strengths in high redshift giant clumps, but numerical 
simulations of magnetized turbulence suggest that a turbulent 
dynamo can rapidly amplify an initi ally sub-Alfvenic field to an 
Alfve n Mach number of order unity l lStone. Ostriker fc Gammiel 
If this occurs in giant clumps at high z, the escape speed 
should be reduced similarly for them. 
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I2OIOU iTacconi et alll2013[) . In our toy model we as- 
sume that the SFR is roughly constant throughout 
the clump lifetime. This is consistent with cosmolog- 
ical simulations, where the SFR in clumps docs not 
show a systematic variation with distance from the disc 
centre despite the cont inuous clump migration inward 
(|Mandelker et al.ll2013h . We have replaced in eq. ^ the 
instantaneous gas mass Afg by the total initial clump 
mass Mc . This should be a close overestimate as long as 
the clump is still far from its gas depletion time. When 
estimating the depletion time, we will replace in 
eq. ([7]) by 0.5 Mc, which will make tgfr larger by a factor 
of 2. 

The stars that form inject momentum into the re- 
maining gas at a rate pin, in the radial direction. In 
this section we will not distinguish between the energy- 
driven and momentum-driven routes for producing this 
momentum. Since feedback is generally dominated by 
massive stars, on e can appeal to the "old stars" limit 
(jKrumholz fc Dek cl 2010), where we are concerned with 
timescales longer than the ^ 4 Myr lifetime of a massive 
staiQ. In this limit, the number of massive stars gener- 
ating feedback at a given time is simply proportional to 
the SFR, so we can write 



Pin = V^nM^ , = V'inVL . 



(9) 



The velocity Vl characterizes the momentum carried by 
the stellar radiation field, 



L 



(10) 



where L is the luminosity produced by these stars. As 
discussed in Section 13.11 this quantity is 



Vl = 190 km s" 



(11) 



corresponding to an energy production of cVl = 5.7 x 
10^^ erg g~^. The dimensionlcss parameter is the 
momentum injection factor, it measures the multiplica- 
tive factor by which the actual injected momentum is 
higher than that one would obtain if stellar radiation 
were the only source of momentum, and if every pho- 
ton were absorbed only once before escaping. For conve- 
nience, we sometimes express it in terms of an effective 
trapping factor (see below) 



V'in = 1 + ftr 



ap ■ 



(12) 



2.3 Wind Properties 



If the momentum injected by stars is able to raise ma- 
terial to speeds V^, }t Vc, this material may be driven 
off the clump in a steady wind. Observed outflows from 
clumps suggest that Kv/K is of order a few. Given the 



* The old stars limit almost certain ly applies to giant clump s at 
^ ~ 2, since these have % > 4 Myr jKrumholz fc DekelllioTolV It 
probably applies to giant molecular clouds in the local universe 
as well, since, although these have % ~ 4 Myr, the best obser- 
vational estimates of GMC lifetimes arc ~ 30 Myr l lFukui et alj 
[200I); this is significantly uncertain, however. The old stars limit 
does not apply to smaller-scale s tructures seen in our galaxy (see 
iKrumholz, Dekel fc McKeell2012l for a more thorough discussion) . 



available supply of momentum, eq. ^ , the actual mo- 
mentum of the wind is 



(13) 



where -^oj is the ejection efficiency representing the 
fraction of the injected stellar momentum that goes into 
the wind. The ejection efficiency could be lower than 
unity if some material is raised to speeds below Vc and 
thus does not escape. Similarly, if the momentum injec- 
tion is spatially distributed rather than point-like, there 
may be some cancellation of the momenta injected at 
different positions, again producing -^oj < 1- We define 
for simplicity the overall momentum efficiency factor in 
driving the wind. 



Pw 

L/c'' 



(14) 



as the factor representing the ratio of actual momentum 
in the wind to the momentum carried by the radiation 
when each photon is counted once. This quantity can 
be addressed observationally, as opposed to ■f/^in and "iAoj • 

The wind mass fiow rate can be extracted from the 
wind momentum and velocity via 

MwKv = Pw = ipwVi^M^ . (15) 
The mass loading factor of the wind is 

The corresponding timescale for mass loss by outflow is 
then 



t^ 



Me 



(17) 



For T] of order unity, this timescale is comparable to the 
SFR timescale, eq. ([5]). 

Note from p^ = M^V^ that the mass- loss rate for 
a given momentum budget is maximized if the wind 
velocity factor 



Vc 



(18) 



is as close to unity (from above) as possible, i.e. if the 
ejected gas is raised to the lowest possible velocity con- 
sistent with escape. If i/ is a constant determined by the 
stellar momentum-driven ejection mechanism and is in- 
dependent of the clump escape velocity (as indicated 
observationally, 1/^3, see Sj5]), then 7] is inversely pro- 
portional to Vc , namely the wind mass- loading factor is 
larger for less massive clumps. 

The extr eme limit of this phenom en on, a s con- 
sidered by iFall. Krumholz fc Matzned ()2010[) and 
IKrumholz fc Dekell (|20lOr is an explosive ejection, 
which occurs when the momentum or energy injection 
is sufficient to raise the entire gas mass to speeds ^ Vc 
in a time of order iff. In this case the feedback is likely 
to sweep up all the material in the clump and eject 
it explosively, halting any further star formation. The 
condition for this to occur is that Pw^ff ~ M^Vc, which, 
using eq. (jl5p . reduces to the condition that 
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1 V, 



(19) 



(Mg/tff) ^" 1/l ■ 

As pointed out bv lKrumholz fc Dekell (j2010[ ). this con- 
dition is extremely difficult to satisfy in giant clumps. 
Observations constrain eff ^ 0.01, and, as we will see 
below, Vc/Vi, ^ 0.5 for giant clumps. Thus achieving ex- 
plosive ejection requires either that ^ 10 — 100, that 
eff in giant clumps exceed the observationally-inferred 
values of ^ 0.01 by a factor of ~ 10 — 100, or some 
combination of both. 



2.4 Clump Depletion 

Given our calculated mass loss rates, we can also con- 
sider the implications for the star formation efficiency 
and lifetime of star forming clumps. First consider a 
clump that never experiences explosive ejection, and is 
simply eroded by a combination of star formation and 
a steady wind at constant rates. The star formation ef- 
ficiency at the time of gas depletion will be 



1 



Ah + A/« 



1 



(20) 



Conversely, a clump that has no steady wind ("ipd — 0) 
and also does not satisfy the criterion for explosive dis- 
ruption, eq. (jl9p . will eventually turn itself completely 
into stars, or will undergo explosive disruption at late 
times when the amount of gas is reduced to the point 
where the clump is no longer in the old stars limit. By 
this point, however, it will have already turned the great 
majority of its mass into stars. A more realistic scenario 
is that a clump experiences a steady wind during its life 
and its star formation efficiency is given by eq. (PU)) . 

One can compute clump depletion lifetimes in an 
analogous manner. In the case of depletion by a steady 
wind, the lifetime is 



^dcp 



Mc 



AL ^ 



(21) 



where Af^ = EvtA/c is the final stellar mass, with e^, from 
eq. (|20l) . For the SFR that enters the last equality of 
eq. (HD) we have replaced A^c in eq. ^ by 0.5A'/c, to 
represent the average between the initial gas mass of 
Ale and the final zero gas mass at depletion, and thus 
refer to the characteristic SFR during the period from 
the onset of the wind to depletion. With ij ^ 1, this is 
much larger than tg. In the case of explosive disruption 
the depletion time of the explosive phase is simply of 
order tg. 



2.5 Clump Migration versus Depletion 

In the case of giant clumps in high redshift galax- 
ies, the lifetime may also be limited because after 
some period a clump will migrate into the galactic 
centre following angular momentum and energy loss 
by torques from the perturbed disc, clump-clump in- 
teraction and dynamica l fraction. The time required 
for this to happen is (Ibekel. Sari fc Ceverinol 120091 : 
ICeverino. Dekel fc Bournaudll2010[) 



"-mig 

where 



Vd 



(22) 



(23) 



is the disc crossing time, Rd is the characteristic disc 
radius and Vd is the characteristic disc circular ve- 
locity. The quantity S is the mass fraction in cold 
disc within the disc radius, which at the cosmologi- 
cal steady state is 5 ~ 0.33 . The Toomre parameter 
is Q ~ 0.68 for a thick disc (iGoldreich fc Lvnden-Belll 
ll965HDekel. Sari fc Ceverinol2009ll . The migration time 
is thus comparable to the orbital time at the outer disc. 
If we approximate id — 3iff , assuming that the clumps 
are overd ensities of ^ 10 with r espect to the background 
disc (e.g. ICeverino et al.l [20121 ). we get 



24fff ~ 12eft(l + r;)idop. 



(24) 



The true clump lifetime will be the lesser of tdcp and 
imig, and the corresponding star formation efficiency 
will be the lesser of and 



A/* 

aT 



mig 



ts 



24 eff , 



(25) 



where the last equality assumes again id — 3tff. This 
expression for e,^mig is valid when tmig is significantly 
smaller than tdep, so the approximation Alg ~ Ad^ in 
the SFR is good. If imig and tdep are comparable, then 
a better approximation for e*^niig should be smaller by 
a factor of ^ 2. From eq. (Hj), if eff ~ 0.01, we learn 
is expected to be smaller than idcp as long 
as (1 + 77) < 8.3. When tmig is the shorter timescale, 
namely when rj is of order unity, the clump bound mass 
fraction remaining after outflow mass loss at the end of 
the migration is 



M, 



c jinig 



A^c 



1 - 



AJ f ■ 



Mr. 



1 - ??e^ 



mig ■ 



(26) 



For 77 ~ 1 this is a significant fraction of the original 
clump mass. 

2.6 Clump Mass Growth during Migration 

Accretion onto the clumps, including clump mergers 
and possibly tidal stripping of the clumps, arc signif- 
icant during the clump migration inward, and may ac- 
tually be the dominant effect in the evolution of clump 
mass. As the clump spirals in toward the disc centre it 
accretes matter from the underlying disc. An estimate 
of the accretion rate is provided by the entry rate into 
the tidal (Hill) sphere of the clump in the galaxy, i?T, 

A/4c=^a/9d(^-RT')CTd. (27) 

Here pd is the density in the cold disc (gas or young 
stars), ttRt^ is the cross section for entry into the tidal 
sphere, and (Jd is the velocity dispersion in the disc rep- 
resenting here the relative velocity of the clump with 
respect of the rest of the disc. The parameter a is ex- 
pected to be of order unity and smaller. 

The tidal or Hill radius i?T about the clump is 
where the self-gravity force by the clump balances 
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the tidal force exerted by the total mass distribu- 
tion ill the galaxy along the galactic radial direc- 
tion. If the disc is in marginal Toomre instability with 
Q ~ 1, this is the same as the Toomre radius of the 
proto-clump patch that cont racts to form the clump 
(iDekel. Sari fc Ceverinol[20n9l) . 



i?T ^ 0.5 5Ra , 

where the clump mass is given by 

2 



Me 



i?Ty 
rJ 



(28) 



(29) 



with Md referring to the mass of the cold disc. Also 
when Q ~ 1, the disc half height h is comparable to 
Rt, 

' - (30) 



i?d 
and 

,5- 



(31) 



We can now evaluate the timescale for clump 
growth by accretion, tac, using eq. pT]). We insert 
Rt from eq. ([28|) . write pd = Md/{2TTRd^h), and use 
eq. (pO|) for h to obtain 



Me 
Mar 



2 

-td. 

a 



(32) 



With a ~ 1/3 (see below), the timescale for doubling 
the clump mass by accretion is roughly an orbital time, 
comparable to the migration time. 

The parameter a represents the fraction of the mass 
entering the tidal radius that is actually bound to the 
clump. If the clump collapses from an initial patch of 
radius Rt, the particles enter the tidal radius with a ve- 
locity distribution that is similar to that of the overall 
disc, namely with a standard deviation ad, and a dis- 
tribution of kinetic energies per unit mass about 0.5 cr^ 
in the clump rest frame. Using eq. ((28)) . eq. p9)) and 
eq. dSD), the binding potential of the clump at Rt can 
be crudely estimated by 



(33) 



GMe 

Rt 

so a significant fraction of the particles entering the tidal 
radius are expected to be bound. 

One can estimate a by referring to the particles 
that actually hit the clump, of radius i?e < ^^t, and 
are bound there. This requirement puts an upper limit 
on the particle impact parameter b prior to entering 
Rt such that the focusing of the orbit would bring the 
particle into i?e with a velocity smaller than the escape 
velocity from the clump Vc at Re- Angular-momentum 
conservation yields b ~ RcVc/crd- Using eq. ^ for Vc 
and eq. ((33)) for ad we obtain Vc/ad — {Rt / RcY^^ ■ 
Therefore in eq. (P7)) 

a ~ b^/RT^ ~ Rc/Rt ■ (34) 

With a typical con traction factor of Rt/Rq — 3 
(jCeverino et al.|[2"012[ ). the estimate is a ~ 1/3. 



For a uniform disc, the relevant density pd is the 
mean density in the disc, pd, within a cylinder of ra- 
dius Rd and height 2h, and a is the same through- 
out the disc. However, for an exponential disc with 
an exponential radius i?exp, the local density pd is ~ 
(0.28, 0.46, 0.70) pd at (3, 2, 1) i?oxp, respectively. There- 
fore, if one uses the mean density in eq. (P7)) . the effec- 
tive value of a in the outer disc could be slightly smaller 
than estimated above. 

The clump growth rate is further enhanced by 
mergers of clumps as they spiral in. If the clumps 
conta in 20% of the disc mass ((Dekcl. Sari fc Ceverind 
[2009I) . and if we require binding when the clump centres 
are at a distance of 2i?e from each other, an analogous 
estimate to eq. gives that the timescale for growth 
by mergers is roughly tmor — (5/2)iac- On the other 
hand, tidal stripping may become more pronounced at 
small radii, which may slow down the mass growth 
rate at the late stages of the migration. We expect the 
timescale for mass growth to be comparable to the crude 
estimate in eq. ((32|) . 

The above estimates for <ac and imor are indeed 
consistent with the findings from liydro-cosmological 
simulations, in which outflows by stellar feedback are 
weak by construction. In these simulations, the clump 
mass is found to be roughly i nversely proportional t o 
distance from the disc centre (jMandelker et al.ll2013[ ). 
When following individual clumps as they accrete, strip 
and merge during migration, they indeed grow in mass 
on a timescale that is comparable to the migration 
timescale. An effective value of a ^ 0.33 in eq. ([5^ 
seems to provide a good fit to the overall clump mass 
growth rate in these simulations (with negligible out- 
fiows). 



3 MOMENTUM BUDGET 

Having developed a basic framework for how the prop- 
erties of star-forming clumps - their star formation ef- 
ficiencies, lifetimes, and outflows - depend on the mo- 
mentum injected by stellar feedback, we now turn to 
various possible feedback mechanisms. Our goal is to 
understand the momentum injection efficiency ^pi^ for 
each mechanism. We note t hat a simi l ar bu dgeting exer- 
cise has been carried out bv lMatznerl ()2002[ ) for Galactic 
giant molecular clouds, and our general approach will 
follow his, applied to a quite different context. 

Unless otherwise stated, all the values below are 
derived using a STARBURS t99 ()Leitherer et al.l 119991: 
IVazguez fc Leithereij I2005D calculation for continuous 
star formation, with all parameters set to their default 
values except that we use an IMF upper limit of 120 Mq 
instead of the default of 100 Af©, though the difference 
in most quantities is smaU (< 10%). We use 120 Mq 
because it is the largest mass for which evolutionary 
tracks are available in STARBURSt99, and because ob- 
servations indicate that the IMF extends to at least thi s 
mass if not significantly higher ([Crowther et al.|[2010[ ). 
We evaluate all quantities at a time of 100 Myr after 
the start of star formation, but since luminosity and all 
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other quantities are slowly varying at times > 4 Myr 
(which is the essence of the old stars limit), different 
choices of age in the range 10 — 300 Myr make a dif- 
ference of at most a few tens of percent. Given this 
uncertainty, wc give all results to two significant digits 
only. 



3.1 Stellar Radiation 

Stellar radiation pressure has received a great 
deal of attention recently as a potential mecha- 
nism for disr upting giant clumps, both in ana- 
lytic models (iMurrav. Quataert fc ThompsonI 120101 : 



i;*/(lM0 yr 1 kpc ^), and Sgas.o 



Sgas/(1 g cm 2). 



iKrumholz &: Dekelll2010^ and in numerical simulations 
(iHopkins, Quataert fc Murravl 20111: Hopkins et al 



20121: iHopkins. Quataert fc Murravl 120121: iGenel et al 



2013) ■ As stated above, our STARBURSt99 calculation 



gives Vl = 190 km s ^ as the momentum budget of the 
direct radiation field. 

Recall our definition of tpin = 1-1- /trap as the ratio 
of the momentum actually injected into the gas to that 
which would be imparted by direct radiation pressure 
alone. Thus V'in = 1 corresponds to a flow that receives 
no momentum from any source but photons that are 
absorbed once and then escape. If the radiation emitted 
by stars is trapped by the high optical depths of a dust 
layer, it is possible that the trapped photons will be 
absorbed or re-emitted multiple times before escape. 
If this occurs, radiation energy can build up, and 
the adiabatic expansion of the radiation-dominated 
region can result in a larger momentum transfer to the 
gas. Some analytic and numerical models of radiation 
pressure-driven feedback assume that this effect will 
produce a value of /tra p ~ t, where r is an approximate 
infrared optical depth (|Murrav. Qu ataert fc Thompson 
2010[ IHopkins. Quataert fc Murrav. ■2012i: iGenel et al 
2012t) . while others assume that radiation-driven 
flows are strictly mo mentum-limited, with -^in 
never exceeding a few (IKrumholz Sz Matznei 20091: 



Fall. Krumholz fc Matznerl 120101: IKrumholz fc Dekell 
20101). due to radiation R ayleigh- Taylor instability 
( Jacauet fc Krumholzj [20TT1 ). This instability punches 
holes in the gas that allow photons to leak out, and 
prevent the buildup of adiabatic radiation-dominated 
regions. 

Recent radiation-h ydrodynamic simulations by 
IKrumholz fc ThompsonI ([2012 , 20 13) have significantly 
clarified the matter. They show that, in the case of 
a radiatively-driven wind, the trapping factor obeys 
/trap ~ 0.5r,, where r* is the optical depth evalu- 
ated using the opacity at the dust photosphere, not 
the far higher opacity found deep in the dust gas 
where the ra diation temperature is high e r, as as- 
sumed e.g. bv IHopkins. Quataert fc Murravl (|2012t) . In 
the old stars limit, for an object of SFR per unit 
area E<. and gas surface dens ity Egas this is given by 
(jKrumholz fc Thompsoi]|2013D 



T, = 0.01 {cVl^)li^ K(o Sgas,0 , 



(35) 



The values of S* and Sgas to which we have scaled 
are typical of observed giant clumps, as discussed be- 
low, and thus we typically have <C 1, and there- 
fore /trap.rad ~ 0. We coucludc that the very large 
trapping fact ors of 10 — 50 as sumed in certai n sim- 
ulations (e.g. Qppenheimer fc Dave 2008 ; Gen cl et al.l 
120121: IHopkins. Quataert fc Murravll2012f ) are unrealis- 
tic. 



Note that the energy in the wind is 



1 



1 , 



I Ac 



(36) 



Thus, as long as ■0^ is of order a few and Kv ^ c, 
only a small friction of the photon energy is used to 
drive the outflow. This is indeed momentum-conserving 
rather than energy-conserving driving of the outflow. 



3.2 Photoionized Gas 

In Galactic giant molecular clouds, the pressure 
of photoionized gas is likely the dominant feed- 
back mech anism that limits the star formation effl- 
ciency fe. g. lWhitworthlll979l: IWilhams fc McKed 



Matznerl 120021: IKrumholz. Matzner fc McKed 



1997 



2006 



where (cVi) 



10 



(cVi)/(lOiOi0/(Mo yr-i)), S,,o 



Goldbaum et al.ll2011l) . Photoionization raises gas to a 
nearly flxed temperature of roughly 10^ K, and that 
temperature is maintained by radiative heating and 
cooling processes. As the gas expands, it transfers mo- 
mentum to the surrounding medium. By integrating 
over the IMF and using a similarity solu tion to compute 
the evolution of expanding H ii regions. [Matzneil (j2002|) 
estimates that this mechanism injects momentum at 
a rate Vin — 260 km s"^. While this is probably the 
dominant feedback mechanism for GMCs in local galax- 
ies, it is likely to be unimportant for giant clumps, for 
the simple reason that such clumps have characteristic 
speeds significantly higher than the ionized gas sound 
speed, c\ ~ lOkms^^ (Knmiliolz & Dekcl 2010). As a 
result, ionized gas will be unable to expand and trans- 
fer momentum to the cold gas . Recent numerical simula- 
tions confirm this conjecture (jPale. Ercolano fc BonneUl 
120121) . We may therefore disregard this mechanism for 
giant clumps at high redshift. 



3.3 Protostellar Outflows 

Protostars drive coUimatcd hydromagnetic outflows 
with launch speeds comparable to the escape velocity 
from stellar surfaces, typically ^ lOOkms"^ for proto- 
stars with radii larger than those of main sequence stars. 
The wind material shocks against and mixes with the 
surrounding dense molecular gas. Because the environ- 
ment the winds encounter is very dense, and the shock 
velocity is not high enough to heat material past the 
~ 10^ K peak of the cooling curve, the post-shock gas 
rapidly cools via radiation, so there is no significant adi- 
abatic expansion phase. Starburst99 docs not include 
protostellar winds, so we adopt an estimate of 

yps~40kms"^ (37) 
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from lMatzneil (|2002[ ). 



3.4 Supernovae 

From our STARBURST99 calculation, supernovae occur 
at a rate 



0.012 {AU/Mq yr-i) , 



(38) 



carry an energjQ ~ 5.8 x 10^^ ere g ^, and carry a 
momentum flux that corresponds to 



Vsn ^ 48kms-i . 



(39) 



This represents only a lower limit on the true mo- 
mentum injected by supernovae, because the extremely 
high post-shock temperatures produced when super- 
nova ejecta encounter the ISM guarantee that radiative 
cooling is inefficient at the early stages. As a result, 
supernova remnants experience an energy-conserving 
phase when they are small, then begin to cool radia- 
tively only once adiabatic expansion lowers their in- 
ternal temperatures sufficiently. During the adiabatic 
Sedov- Taylor phase and the pressure-driven snowplow 
phase that follows it (during which the remnant inte- 
rior is partially radiative), the radial momentum carried 
by the swept-up material increases. T his process has 
been studied by numerous authors (e.g. Chevalieil 1974 ■ 
McKee & Ostriker"l977l; ICioffi. McKee fc Bertschingeij 
1988; Thornto n et al. 1993) j for a uniform medium 
Thornton et al.l find that the asymptotic momentum of 
a supernova remnant is 



1.7 X IQ-'^E, 



,43 7pl3/14 ^-0.25 



51 



g cm s ^ , (40) 

where E^i is the energy of the supernova in units of 10^^ 
erg and ni is the ambient number density in units of 10 
cm""^. A simple estimate for the supernova momentum 
budget including the adiabatic phase is simply obtained 
from psn, which is given by Psn from eq. (|40)) multiplied 
by the supernova rate from eq. ([55)) . This gives 



sn.adiab 



1100nj;°-2^ km s" 



(41) 

assuming the energy of a single supernova is E^i = 1 (as 
assumed in the STARBURSt99 calculation as well). Since 
the momentum input from a single supernova remnant 
is very close to linear in E^i, the total momentum bud- 
get is not significantly affected by the manner in which 
the supernovae are clustered. 

However, we caution that this calculation is for 
a simple, one-dimensional uniform medium. As the 



^ Note that V and here are the energy and momentum per 
unit mass of stars formed, not per unit mass of stars that actually 
end their lives as SNe. 

^ The publicly-available version of STARBURST99 does not calcu- 
late the supernova momentum flux. We have modified it to do 
so, using the same assumptions STARBURST99 adopts in order to 
compute the supernova energy injection rate and mass return, 
i.e. all stars with an initial mass above 8 Mq end their lives as 
supernovae with identical energies of 10^^ erg. They all leave be- 
hind as remnants 1.4 Mq neutron stars, so the mass of the ejecta 
is simply the final stellar mass (smaller than the initial mass due 
to wind losses) minus the remnant neutron star mass. 



iKrumholz fc ThompsonI (|2012l . l2013f ) results for radia- 
tion pressure show, this assumption can be deeply mis- 
leading about how effectively energy is converted into 
momentum in a real three-dimensional medium where 
instabilities can occur. It is therefore best to regard 
eq. (|4ip as representing an upper limit. Determining 
where reality lies in between this value and the lower 
limit represented by eq. (j39p requires numerical simu- 
lations capable of following instabilities into the non- 
linear phase. Although suc h simulations have begun 
to ap pear in the literature (jCreasev. Theuns fc Bowerl 
I2013I) . the problem remains far from fully solved. 



3.5 Main Sequence and Post-Main Sequence 
Stellar Winds 

The winds of main sequence and post-main sequence 
stars carry an energy and momentum content ~ 
1.5 X 10^^ erg g~^ and 



Vr 



ms.dir 



140 kms 



(42) 



respectively. While they therefore carry slightly less mo- 
mentum than the stellar radiation field, at least some of 
the winds are launched at velocities large enough that 
the post-shock gas may have long cooling times. (This 
is in contrast to the much slower protostellar outflows.) 
As a result, it is plausible that stellar winds could expe- 
rience an adiabatic phase like supernovae, and enhance 
their momentum transfer that way. In this case the mo- 
mentum provided by winds will have an additional term 
that we can write as fudVh- 

The main idea of the classical stella r wind bub- 
ble model of ICastor. McCrav fc Weaved ()1975D and 
[Weaver et"all (|1977[ ) is that f.^iVr, » T4...H;r. On the 
other hand, this mechanism will not operate if wind gas 
is able escape from a star-forming clump without en- 
training significant mass, or if it undergoes rapid cooling 
by mixing with cooler, dense gas that brings its temper- 
ature low enough for radiative losses to become rapid. 
Observations of a few nearby H ii regions have been 
able to address this question directly by using x-ray 
observations to p robe the energy density of the shock- 



heated gas. Both iHarper-Clark fc Murravl ( 2009 ). who 
study the Carina Nebula. and lLopez et al.l ( 2011 ) , who 
study 30 Doradus, find that the luminosity of the x-ray 
emitting gas implies that the pressure exerted by this 
gas is weaker than that exerted by photoionizcd g as, a 
result highly inconsistent with an energy-driven flowQ 



^ Contrary to the resu lts of iLopez et al. I l l20lj) , 

iPelleerini. Baldwin fc Ferlan argue that the x-ray 

emitting gas pressure i n 30 Doradus is actually high er than 
the ionized gas pressure. iPellegrini. Baldwin fc Ferlandf s results 
differ because they assume that the x-ray emitting gas is confined 
to a small volume. Since the x-ray luminosity is proportional to 
the emission measure of the emitting gas, which is the integral of 
the square of the electron density along the line of sight, if one 
assumes that the line of sight length is much smaller than the 
transverse size of the region being observed, the density and thus 
pressure that one infers for a given observed luminosity rises 
proportionately. For our purposes, however, this distinction is 
irrelevant. If one assumes that the hot gas is confined to a small 
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For our fiducial estimate in this paper we adopt 
Lopez et al.'s measured mean value /ad = 0.3, so that 
the net amount of momentum ejected by the winds from 
main-sequence stars is 



230 + 1350 n^'"-^^ km s 



^ms — Kns,dir + fadVh 



200 kms" 



(43) 



However, we caution that none of the observed regions 
have conditions close to those of high rcdshift giant 
clumps. While stellar wind gas is momentum- and not 
energy-driven in the local universe, a giant clump could 
be considerably harder for hot x-ray gas to escape. It is 
therefore conceivable that, under the conditions found 
in high redshift giant clumps, stellar winds represent an 
energy-driven feedback. In this case it is likely that stel- 
lar wind and supernova bubbles would simply add to- 
gether to produce an adiabatic shell driven by the com- 
bined effects of both. The result is to increase the adi- 
abatic energy budget by roughly 25% compared to su- 
pernovae alone. Adopting this simple estimate, we find 
that if winds are adiabatic then, in conjunction with su- 
pernovae, the net momentum contribution of the winds 



ms.adiab 



275nj;°-25 km s"\ 



(44) 

Note that this estimate implicitly assumes that the tem- 
perature and cooling are determined by the significantly 
larger energy associated with the supernovae, so that 
stellar winds simply pump more energy into the adia- 
batic bubble without significantly affecting how it cools. 



3.6 Total Momentum Budget 

Combining all the mechanisms we have enumerated (ex- 
cluding photoionized gas for the reasons stated above), 
we see that the momentum budget for the case of a 
purely momentum-driven outflow is expected to be 



480 km s"^ (45) 
The corresponding trapping factor is rather small, 



V'in = 1 + ./tr 



ap 



2.5. (46) 

By "purely momentum- driven outflow" we refer to the 
case where radial momentum is conserved and nei- 
ther supernovae nor stellar winds experience an energy- 
conserving phase during which their radial momentum 
is significantly boosted. This can be interpreted as a 
lower limit for from stellar feedback. 

If we assume that there is a significant adiabatic 
phase for supernovae and main-sequence winds, and 
that the resulting momentum injection is near the up- 
per limit derived in the uniform medium case, then we 
obtain an upper limit for the net momentum budget of 



in.adiab 



^■ad + Vps + Kn,adiab + Vi 



ms.adiab 



fraction fx of the observed volume, the pressure Px inferred 

— 1/2 

from a given luminosity varies as Px ' , but the energy 

1 /2 

content of the hot gas, which v aries 

fxPx, falls as . Thus if 
iPelleerini. Baldwin &: Ferlandl 's conjecture about the geometry 
is correct, that implies even more strongly that feedback from 
stellar winds cannot be energy-driven. 



1.2 + 7.177.7 



(47) 
(48) 

For massive clumps of tff ~ 7Myr this is ipia — 6.4. 

The actual contribution of adiabatic supernova 
feedback and the corresponding value of ipin between 
the above lower and upper limits is a matter of a 
saturated state of fully nonlinear instabilities, which 
should be determined by appropriate numerical sim- 
ulations. The two most relevant publications to date 
on t his subject are Hopkins. Quataert fc Murray! (|2012l) 



and ICreasev. Theuns fc BoweiT i 2013f ). For the former, 
if we examine the runs without the subgrid radiation 
model where supernova feedback dominates, the mass 
loading factors are '--^ 1 — 5 and the wind terminal 
speeds are a few hundred km s~-^, implying ip^ values 
of a few; combined these sugges t Vin ~ 5. Similarly, 
ICreasev. Theuns fc Boweil ()2013[ ) report a mass load- 
ing factor r] and a wind thermalization parameter tjt 
(which measures the fraction of supernova energy that 
goes into outflow; see their equation 5), both as a func- 
tion of galaxy properties. With some algebra, one can 
show that = 5.3(7777t)^^^, and using their fitting for- 
mulae for galaxy surface densities ~ 10— 100 Mq pc~^, 
appropriate to giant clump galaxies, gives ip^/ ~ 1 — 3 
for supernov ae alone. In summary, the numerica l re- 
sults of both [ Hopkins. Quatae rt fc Murray! ()2012l) and 
ICreasev. Theu ns & Bower (2013) suggest that Vin is 
likely to be roughly halfway between our upper and 
lower limits. 

A potential way to make the supernova feedback 
more effective is by having the gas density in the super- 
nova vicinity much lower than the unperturbed density 
within the clump of 7ii ^ 1. A value of ni ~ 10^^ (or 
10~^) in eq. (|48|) would provide a maximum value of 
V'in — 24 (or 41 respectively). The question is whether 
the other types of momentum-conserving stellar feed- 
back could generate such a low-density regime prior to 
the supernova explosion. We keep this mechanism out- 
side the scope of the present paper. 



4 IMPLICATIONS FOR HIGH-REDSHIFT 
GIANT CLUMPS 

4.1 Star Formation and Outflows 

In the preceding two sections, we developed a general 
framework to consider the evolution of clumps as they 
migrate, accrete, form stars, and lose mass due to star 
formation feedback, and we derived an estimate for the 
momentum budget of the feedback that drives clump 
winds. We now combine these results to draw conclu- 
sions about the typical evolutionary path taken by giant 
clumps. 

The structure and dynamics of a clump are char- 
acterized by two quantities, e.g., its characteristic ve- 
locity and its free-fall time, Vc = 100kms~^V2, and 
tff = lOMyrtffio, as defined in eq. 

* The relations to the clump mass and radius are V2 — 
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Following the earlier discussion, the physics of out- 
flow from giant clumps can be characterized by three 
dimensionless parameters, e.g., eg, v and '0w = V'ojV'm- 
These are the SFR efficiency eff = O.Oleff _2, the wind 
velocity with respect to the clump escape velocity, 
u = V-^fi/Vc = Si's, and the wind momentum with re- 
spect to the radiation momentum, ^pv, = Pw/iL/c). For 
the values of eg and ipin we have theoretical predictions. 
For the values of ly and tpej , unity is a lower and an up- 
per limit respectively, but we do not have a theoretical 
prediction concerning how much they actually deviate 
from unity. Motivated by observations (see below), we 
assume that these deviations are by a multiplicative fac- 
tor of order one or a few. We define i/j^ = 2.5V'w,2.5- 
The reference value of ipw — 2.5 may refer to the case 
of pure momentum-driven outflow -^in ~ 2.5 and maxi- 
mum ejection of tpcj ~ 1, or to a case including adiabatic 
supernova and stellar-wind feedback with ipi^ — 5 but 
with some losses in the ejection, ipcj — 0.5. The maxi- 
mum value of when adiabatic supernova feedback is 
at its maximum and the ejection is efficient, is expected 
to be ipv/ ~ 5. 

The SFR and wind mass flow rate are 



Mw ^ 3.2eff,_2Vw,2.5i^3"^^2^ A^oyr"^ ■ 
The corresponding timcscales are 



10 ■ 



(49) 
(50) 

(51) 
(52) 



The mass loading factor, eq. ([T6|) . is 

V - V'w,2.5^w;400 = l-33'0w,2.5'^3"^1^2"^ > (53) 

where = 400 km s~^ 1/^.400 ■ With the fiducial val- 
ues adopted here for momentum-driven stellar feed- 
back from typical clumps one expects steady winds 
with mass-loading factors of order unity. The maximum 
value, when adiabatic supernova feedback is included, 
is expected to be ?7 ~ 2 — 3. 

If the clump were allowed to deplete all its gas, the 
final star formation efficiency would have been 

e. = (l + 77)-i = 0.5(1 + 77)^1 (54) 

at the clump depletion time of 

(55) 

where (1 + r/)2 = (1 + ri)/2. 

The approximate values for and tjop are valid 
when e^, deviates significantly from unity, and where 
clump disruption is by gradual erosion rather than sud- 
den explosive destruction. The criterion for explosive 
disruption, eq. p^ . is simply 



tdcp = 1 Gyr (1 + 77)2 ^Eff _Vff ,10 1 



LlSMg/^Rj and tw ^ 
1.27M9.5y„"^ where M9.1 



A/c/1O«-5M0 and _Ri 



-Rc/lkpc. Also rii 
0.21M9.5/?j;^gcm 



2.2t,Q . The surface density is S 



eff,-2'/'w,2.5^2'' ^ 20. (56) 

This is similar to eq. 9 of iKrumholz fc Dekell (|2010[ ). 
where the considerations were qualitatively similar 
though not exactly the same numerically. As noted 
there, if the clump is a typical Toomre clump with 
Vc ^ lOOkms"^, explosive disruption occurs only if ei- 
ther eff or ipv, are significantly larger than their fidu- 
cial values, namely either the SFR is much more ef- 
ficient than implied by the local Kennicutt relation, 
or the momentum-driven feedback is much more effi- 
cient than available in the momentum budget evaluated 
above. Otherwise, with the adopted fiducial values for 
these quantities, explosive disruption does not occur, 
thus validating the steady-wind approximation used. 



4.2 Clump Migration 

In comparison, the clump migration time is 



8id ^ 260Myri?d,7V'd 



-1 
,200 ' 



(57) 



where the disc is characterized by — Rd/Vd with 
Rdj = i?.d/7kpc and Vd,2oo = Ki/200kms-i. The fidu- 
cial values of Rd and Vd are d e duced from observations 
at z - 2 (|Genzel et al.ll2006l [20081 ). but the disc dy- 
namical time can also be derived from the virial radius 
and velocity using the virial relation and the spherical 
collapse model and assuming a constant spin parame- 
ter A for halos and conservation of angular momentum 
during gas collapse within the dark matter halo. 



td 



0.07 Ao.07 0.15 tHu 



bble : 



(58) 



which at z = 2, where ^Hubble — 3.25 Gyr, gives td ~ 
33 Ao.07 Myr. 

The relation between the depletion time and migra- 
tion time is 



tmig/tdcp ^ 0.25 (1 + lf)2 tS-2 , 



(59) 



where we have assumed td — 3tff for the dynamical 
timescales in the disc and in the clumps. If tniig ^ tdcp, 
The maximum star formation efficiency possible before 
the clump reaches the galactic centre is 

e,,„ig ~ 0.24eff^_2 • (60) 

The clump bound mass fraction remaining at the end 
of the migration is 



c,mig ^ -^^w^mig 



Mr. 



Mr. 



1 - 0.2477eff 



(61) 



^ and 1 e cm 



: 4800A/, 



Pc 



These are for tmig significantly smaller than tdcp, namely 
for 77 of order unity. The estimate for e^^mig is an overes- 
timate by up to a factor of order 2 because we assumed 
here Mg ~ Mr- For the same reason, the expression for 
Mc inig/Mc is an underestimate. With the fiducial val- 
ues adopted here, the clump reaches the centre while 
still holding on to a significant fraction of its original 
mass, and most likely still gas rich. 



4.3 Clump Mass Evolution 
Eq. (I32|), with td ~ 3tff, yields 
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iac ^0.18Gyrao;j3iff,io, (62) 

where 0:0.33 = a/0.33. Thus, with a = 0.33, the 
timescale for doubhng the clump mass by accretion is 
~ 6td ~ 0.18 Gyr, which is shghtly smaller than the 



migration time, t„ 



8td 



0.24 Gyr. 



With the fiducial values for momentum-driven 
winds, the mass growth rate as estimated in eq. (|62p 
is faster than the outflow rate and the SFR, eq. ([5^ 
and eq. (|5ip . implying that the accretion more than 
compensates for the mass loss by outflows, making the 
clumps actually grow in mass as they migrate inwards. 
This implies in particular that the adopted estimate of 
migration time remains a good approximation and may 
even be an overestimate. 

The evolution of clump mass M{t) under accretion 
and outflows, starting from an original mass Mc at t — 
0, is governed by 

M = Mac - Mw . (63) 

What makes the integration of this equation simple is 
that the two terms on the right hand side both scale 
with M/tff. First, 



Ma 



a tff M 
2 <d tff 



(64) 



where a and iff /id are approximated as constants, the 
latter being determined by the clump collapse factor 
from the original protoclump patch in the disc. Second, 



M^ 



f 



(65) 



where /g is the star-forming gas fraction in the clump, 
approximated as constant. Integrating we obtain 

M(0 = Mce^*/*«, 7 = 0.5a(tftAd)-??eff/g. (66) 

With our fiducial values (a = 0.33, tff/td = 1/3, 
r, = 1, eff = 0.01, /g = 1) we have 7 = 0.045. With 
^mig = 8id the growth factor during migration becomes 
M(tniig)/Mc ~ 2.9. It requires a very strong wind of 
T] ~ 5.5 for the mass loss to balance the accretion and 
leave the clump with a constant mass till depletion, 
which in this case may occur before the clump completes 
its migration. For a significant mass loss in a migration 
time, tmig ~ 24tff, 7 in eq. has to be significantly 
smaller than —1/24. With the fiducial value of a = 0.33, 
this requires that rjeff would be larger than its fiducial 
value of 0.01 by an order of magnitude. Alternatively, 
7 could obtain such negative values if the effective a is 
negative, e.g., representing a case where mass loss by 
tidal stripping overwhelms the mass gain by accretion. 
However, the reported significant clump growth in the 
simulations, where both accretion and tidal stripping 
are at play, indicates that the effective a is positive and 
close to the assumed fiducial value. We conclude that a 
net mass loss in the clumps is very unlikely. 



4.4 Other Implications 

The predictions listed above have a few immediate and 
interesting implications. If winds are relatively efficient. 



i.e. ipej ~ 1, then when all types of stellar feedback are 
taken into account one expects giant clumps to experi- 
ence fairly significant steady winds. Purely momentum- 
driven feedback is expected to provide a mass loading 
factor r/ of order unity, and adiabatic supernova feed- 
back can boost it to of a few. We emphasize that the 
significant outfiows hold even though the radiative trap- 
ping is negligible, and even though the clumps do not 
experience explosive disruption in a dynamical time. 

If the clumps were allowed to reach depletion, the 
depletion time would have been on the order of a signif- 
icant fraction of 1 Gyr. However, the clumps are likely 
to complete their migration inward at a shorter time. 
During this migration, the clumps accrete mass from 
the disc and merge with other clumps, roughly doubling 
their mass in one orbital time. 

The fact that the timescale for migration is 
typically shorter than the timescales for star for- 
mation and depletion indicates that the clumps 
complete their migration while still gas rich, 
thus taking par t in the overall "wet" inflow 
within the disc (iForbes. Krumholz &: BurkertI 120111 : 



ICacciato. Dekel fc Genej 120121: iDekel et al.l \20m . 

This "wet" inflow has interesting implications, e.g., it 
naturally leads to a compact bulge (Dekel & Burk- 
ert, in prep.) and could fe ed the central black hole 
(jBournaud et al.ll201lL [2012f ). 

A question often raised is whether the outflows from 
clumps can be the driver of turbulence in the disc, 
the mechanism that maintains the Toomre instability 
at Q ~ 1. A necessary condition is that the power in 
the outflows is comparable to the turbulence dissipative 
loss. The outflow power from Nc clumps is 



A^cMwK, 



(67) 



with Mw and as predicted above. The turbulence 
is expected to decay on a disc dynamical time, so the 
dissipation rate is 



(68) 



One can use from our analysis above Vw = i^Vc, M„ = 
rjMi^, Mi, = effMg/tff. For a clump contraction factor 
c, the dynamical times are related as td = ctff. For 
a Q ~ 1 disc and clumps one can estimate that the 
internal clump velocity and the external disc velocit y 
dispersion are comparable (e.g. ICeverino et al.ll2012[ ): 
they are related by 



V,yal^{n/2)c. 



(69) 



Then the ratio of the rate of energy injection by clump 
winds to energy loss due to decay of turbulence becomes 



-E-dis 



6iVc,5?7eff-2 i^: 



3 "-3 



(70) 



This seems to indicate that there is enough energy 
in the outflows to continuously stir up the disc. 
However, it is likely that a large fraction of the 
outflow energy will be ejected along the descending 
density gradient perpendicular to the disc and not 
injected into the inter-clump medium in the disc 
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plane, thus making the contribution of outflows to 
the disc turbulence only secondary. The gravitational 
gain by the VDI-driven inflow along the potential 
gradient within the disc is a more likely source of 
energy for ma i ntainin g the disc turbulence (e.g. 
Bournaud et al.l l201lt iForbes. Krumholz fc Burkera 



201lHCacciato. Dekel k Genelll2012HDekel et al.ll20m 



5 COMPARISON TO OBSERVED CLUMPS 
5.1 Observed Clumps 

Table [T] lists pioneering estimates of the properties of 
seven giant clumps as observed in five z ~ 2.2 star- 
forming disc galaxies (SFG) using AO spectroscopy fo- 
cusing on Ha at the ESQ VLT as p art of the SINS 
survey (jForster Schreiber et al.l 20091). These data are 
based on Table 2 of lGenzel et al.l (j201lh . with slight re- 
visi ons for the mas s ive clu mps in ZC406690 from Table 

3 of lNewman et al.l (j2012|) . The five galaxies are selected 
to be massive discs of rotation velocities ^ 250 km s^^, 
dynamical masses of more than IO^'^^Mq within the in- 
ner lOkpc, and SFR~ 120 - 290 yr'^. They sample 
the upper end of the SFG population, and therefore the 
most massive giant clumps. 

The galaxies BX482 and ZC406690 are large 
clumpy rotating discs with a prominent 5 kpc ring 
of clumps and star formation. D3al5504 is a large ro- 
tating disc with a central AGN. ZC782941 is a more 
compact rotating disc, showing an asymmetry due to a 
compact clump off the main body of the galaxy, poten- 
tially a minor merger. BX599 is a compact system with 
a high velocity dispersion and a small 3 kpc rotating 
disc. 

The most prominent clumps were identified from at 
least two different maps of Ha velocity channels. Clump 
1 in Table [T] is the dominant clump (A) in BX482, part 
of a ~ 5 kpc ring that includes 3 additional smaller 
clumps. Clump 2 is an average over the 6 off-centre 
clumps (A-F) in D3al5504, none of which is particularly 
dominant over the others. Clump 3, from ZC782941, 
is at the largest distance from the disc centre and the 
brightest in Hq, , while this galaxy shows four additional 
clumps closer to the centre. Galaxy ZC406690 shows 

4 clur nps in a ^ 5 kpc rin g, of which three were stud- 
ied in iGenzel et al.l (|2011| ) and listed here. Clump 4 is 
ZC40669G-C, the brightest in I-band ACS, faint in H„, 
and shows an elongated shape. Clump 5 is ZC406690- 
A, the brightest in Hq, and rather round, compact and 
isolated. Its SFR is high, its stellar population is young, 
and it is gas rich. Clump 6 is ZC406690-B, the second 
in Hq, brightness and rather faint in I-band ACS. Its 
stellar population is rather old, and it is relatively gas 
poor. Clump 7 is an exception, the whole centre of the 
compact galaxy BX599, namely a compact star- forming 
bulge. 

The first group of rows in Table [T] refer to the clump 
structural properties. The second group of rows are the 
observed SFR and wind properties. The third group is 
quantities deduced from the observed quantities. The 



quantities marked by "*" are directly deduced from the 
observations. 

As descr ibed in iGenzel et al.l ()2011h and 

iNewman et al.l (|2012t ). the quoted quantities are 
highly uncertain. They are limited by resolution and by 
modeling assumption s. For example, the fo rmal errors 
quoted in Tabic 3 of E ewman et al l (|2012D are about 
100% for some of the quantities characterizing the 
winds. These pioneering observations should therefore 
serve as preliminary indications only. 

The intrinsic clump radius Rc was determined from 
the HWHM of a Gaussian fit to the appropriate veloc- 
ity channel after subtracting in quadrature the HWHM 
of the instrumental resolution. Since the latter is typi- 
cally 2 kpc, larger than the intrinsic clump radius, the 
estimated i?c is rather uncertain. 

The clump characteristic velocity Vc is derived here 
from the kinematic measurements of velocity disper- 
sion (T and rotation V,ot assuming Jeans equilibrium: 
Vc^ = -I- CCT^). The steep clump density profiles 

dictate c ~ 3 .4 (R. Genzel, pri vate communication), 
and ^ ~ 1.17 (jGenzel et al.ll201lt) . Then the dynamical 
clump mass is derived from = Q-^Vc'^Rc. Thi s gives 
larger masses than derived in lGenzel et al.l (j201ll ) using 
c = 2, the value approp r iate f or an isotropic isother- 
mal sphere. IGenzel et al.l (|201l[ ) evaluated the clumps' 
gas mass from the measured SFR using an adopted 
version of the Kennicutt- Schmidt (KS) law. With the 
recent calibration a t z 2 using CO measurements 
('Ta cconi et al.|[20T3l ) , the Kennicutt relation is not very 
different than eq. ([7]) with eff ~ 0.01, and the estimated 
gas mass using the KS law with the recent calibration 
is similar to the dynamical mass as derived here. Note, 
however that the dynamical mass could be underesti- 
mated if the clump deviates from equilibrium due to 
strong infiows or outfiows. 

The SFR is derived from the Hq luminosity cor- 
rected for extinction, with an uncertainty of about 30%. 
The clump outfiow velocity V„ is estimated from the 
maximum blue shift and width of the broad emission 
compone nt - this i s the main pioneering discovery of 
fCenzel ct~all (|201lD . Its error is about 33%. However, 
by adopting the maximum wind velocity as the charac- 
teristic wind velocity one may overestimate some of 
the calculated wind properties. The mass outflow rate 
Mw is taken as the average of two different crude esti- 
mates using two photodi ssociation case-B rn odels, de- 
scribed in appendix B of lGenzel et al.l (|201lD . Because 
of the elaborate modeling involved, and the different re- 
sults obtained from the different models, this quantity 
is naturally highly uncertain, with an error on the order 
of 100%. This is what makes the current results indica- 
tive only, not to be taken too strictly on a case by case 
basis. 



5.2 Comparison of theory to observations 

The first four clumps, addressed as "typical" clumps, 
seem to be consistent with the fiducial case discussed 
above for stellar momentum-driven outflows. The SFR 



© 2013 RAS, MNRAS OOO.nTOI 



Outflows in Giant Clumps 13 



Table 1. Observed properties of giant clumps form lGenzel et aP l|201ll ). Quantities marked by "*" are deduced relatively directly from 
the observations, while the other quantities are computed from them. 



Clump no. 


1 


2 


3 


4 


5 


6 


7 


Clump name 


BX482-A 


D3al5504-A-F 


ZC782941-A 


ZC406690-C 


ZC406690-A 


ZC406690-B 


BX599 


z* 


2.3 


2.4 


2.2 


2.2 


2.2 


2.2 


2.3 


Rc [kpc]*^ 


1.0 


1.0 


0.8 


1.2 


0.8 


1.2 


1.5 


Vc [km s-i]*^ 


125 


111 


195 


159 


163 


187 


152 


Mc [lO^Mo] J 


3.6 


2.8 


6.9 


6.9 


4.8 


9.5 


7.8 


% [Myr]f 


7.7 


8.7 


3.9 


7.2 


4.7 


6.2 


9.5 


[Mgyr-i]*^ 


12 


3.3 


17 


14 


40 


11 


66 


Vw [km s-i]*_/ 


350 


400 


420 


355 


440 


810 


1000 


Mw [Moyr-lpS. 


12 


3.6 


34 


13 


117 


78 


185 


h 


2.5 


1.0 


1.0 


1.5 


3.8 


0.7 


7.9 




2.8 


3.6 


2.2 


2.2 


2.7 


4.3 


6.6 


vl 


1.0 


1.1 


2.0 


0.9 


2.9 


7.1 


2.8 




2 


3 


4 


2 


6 


34 


14 


Sgaa [g Cm-2] J_ 


0.26 


0.20 


0.77 


0.34 


0.54 


0.48 


0.25 


tdcp [Myr]:;^ 


302 


817 


274 


517 


62 


218 


63 




0.50 


0.48 


0.33 


0.52 


0.25 


0.12 


0.26 


^mig/^dep ^ 


0.63 


0.27 


0.36 


0.35 


1.9 


0.71 


3.7 


^* ,niig ^ 


~ 0.50 


0.24 


0.23 


0.35 




~ 0.12 




A^c,mig/A/cl 


~ 0.50 


0.73 


0.54 


0.68 




~ 0.12 





-Rc = ^HWHM after beam smearing (of HWHM~ 2 kpc) is subtracted in quadrature 
^ Vc'^ = P (V^^^ + ca^), /3 = 1.17, c = 3.4, assuming Jeans equilibrium 
c Mc = G-^Vc^Ra 

iff = Rc/Vc 

^ Mi, = L(Ho)/2.1 X 10**^ crgs^^ extinction corrected 

Vw = (Vbroad) " 2(Tbi.oad 

9 Average of two photodissociation caso-B models. Clumps 5 and 6 arc from lNewman et al.l 11201! 

ea^_2 =A.2Vc-2^AU,io (eq.HSll 
\ V = Vw/Vc 
j r] = M^/Mi, 

Vw = M^V^/{L/c) 
'■ Egas = Mg/(7riJc2), Mg = eff-lM^iff 

idcp = 1000Myr(l+7?)-leff-_l2tfi_io (cq.[53 

= (1 + r])~^ (eg. 1541 relevant when fniig/*dcp > 1) 
° tmigAdep = 0.25(1 + r?)2eff^_2 (cq.® 
^ s*,mig = 0.24eff _2 (eg. 1601 relevant when tmigAdep < 1) 
1 Mc^rnig/ = 1 — '?e*,mig (eq. ET] ignoring accretion) 



efficiency eg _2 is of order unity, the wind velocity is 
2-4 times the clump velocity, ~ 3, the mass loading 
factor J] is about unity and in one case ^ 2, and the 
momentum injection-ejection factor Vw is 2-3 and in 
one case ^ 4, as predicted by the theoretical momentum 
budget discussed above. 

The 5th clump, ZC406690-A, is unusual in terms of 
its high SFR of 40MQyr~-'^, but is still marginally con- 
sistent with the fiducial momentum-driven wind case. It 
shows a marginally high SFR efficiency of eff^_2 — 3.8. 
Its outflow is on the strong side, with 77 ~ 2.9 and 
ijj^ ~ 6; but this is still marginally consistent with the 
fiducial case. However, it is different in the sense that 
its high SFR and low iff yield a short depletion time of 
~ 60 Myr, which is about half the migration time. This 
clump will complete its migration intact only because 



the mass gain by accretion from the disc is expected to 
be larger than the mass loss by outflow. 

The last two clumps seem to be extreme cases of 
strong outflows that are inconsistent with stellar feed- 
back, even when the the adiabatic supernova feedback 
is at its maximum and the ejection into the wind is ef- 
ficient. Clump 6 is the most extreme case in terms of 
outflow. It is the most massive clump, ~ IO^^Mq, its 
SFR efhciency is rather typical, ^11 Mq yr~^, with the 
stellar population rather old, but its outflow is exces- 
sive, with 77 ~ 7 and ip^j — 34. Based on our estimates 
of the momentum budget, such an outflow cannot be 
driven by stellar feedback. Either it requires another 
driving mechanism, the observational estimates are se- 
vere overestimates, or the SFR we measure today is 
substantially smaller than it was when the outflow was 
launched. 
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The bulge clump 7 is an exception, representing a 
whole galaxy rather than a Toomre clump embedded in 
a disc. It has the highest SFR~ 66MQyr~^. It shows 
an outflow with a moderately large rj ~ 2.8 but with 
a very high momentum injection efficiency of '0w — 14. 
As a result, its depletion time of idep — 63Myr is only 
a quarter of its migration time. According to our mo- 
mentum budget, such a high value of ip^j is more than 
what stellar feedback can offer in clumps; again, one 
could avoid this problem if the SFR we measure today 
is lower than it was when the bulk of the outflowing 
material was launched. 

[Andrews fc ThompsonI ()2010D proposed a scenario 
where multiple scattering is possible when the gas sur- 
face density is above a threshold of Sgas-^ ~ 1-1 gcm~^. 
If this was true, and if the gas surface density in the ex- 
treme clumps was sufficiently high, this could have pro- 
vided a possible explanation for the extreme clumps. 
Howev er, the simulations of iKrumholz fc ThompsonI 
(|2013f ) show that multiple scattering does not occur 
even at a very high surface density. Furthermore, we 
note that the gas surface density in all the observed 
clumps is in the range (0.2 — 0.8) gcm~^, below the sug- 
gested threshold value given that the metallicity is com- 
parable to and slightly lower than solar. There seems to 
be a marginal correlation between Sgas and V'w, but it 
does not help us explain the extreme outflows in clumps 
6 and 7. 



6 CONCLUSION 

We have analyzed the outflows expected from star- 
forming giant clumps in high-z disc galaxies that un- 
dergo violent disc instability (VDI). We evaluated the 
outflow properties based on the momentum budget, 
namely the efficiency of momentum injection into the 
ISM per unit star formation rate by a variety of stel- 
lar momentum and energy sources. We then estimated 
the lifetime of the clumps given their VDI-driven migra- 
tion toward the disc centre and the associated growth 
of clump mass by accretion from the disc. 

Our results can be summarized as follows: 

• Most of the mass loss is expected to occur through a 
steady wind over many tens of free-fall times, or several 
hundred Myr, rather than by an explosive disruption in 
one or a few free-fall times, less than ^ 100 Myr. 

Radiation-hydrodynam ics simulations 



bv IKrumholz fc ThompsonI (j2012D and 



IKrumholz fc ThompsonI (|2013D provide a key in 



put to the momentum budget, that radiation trapping 
is negligible because it destabilizes the wind. This 
means that each photon can contribute to the wind 
momentum only once, and the radiative force is limited 
to about L/c. This calls into question other recent 
works that assume a very large trapping factor without 
self-consistently computing it. 

• All the direct sources of momentum taken together 
inject momentum into the ISM at a rate of about 
2.5 L/c. This includes radiation pressure, protostellar 



winds, main-sequence winds and direct injection of mo- 
mentum from supernovae. 

• The early adiabatic phases in expanding supernova- 
driven shells and main-sequence winds, if they operate 
at maximum efficiency, bring it up to a total force of 
5 L/c for typical gas densities in the clumps. An un- 
known fraction of this force is actually used to drive 
the wind, so this can serve as an upper limit. 

• The resulting outflow mass-loading factor is of order 
unity. If the clumps were allowed to deplete their gas 
into stars and outflows standing alone, the depletion 
timcscalc would have been a few disc orbital times, a 
significant fraction of a Gyr, ending with about half the 
original clump mass in stars. 

• However, the clump migration time to the disc centre 
due to the VDI is on the order of an orbital time, about 
250 Myr, so the typical clumps are expected to complete 
their migration prior to depletion. 

• Furthermore, based on analytic estimates and simula- 
tions, the clumps are expected to double their mass in 
a disc orbital time by accretion from the disc and merg- 
ers with other clumps, which overwhelm the mass loss 
by tidal stripping. This high rate of gravitational mass 
growth implies a net growth of clump mass in time and 
with decreasing radius despite the continuous massive 
outflows. 

• From the six disc clumps observed so far, five are con- 
sistent with the predictions for stellar-driven outflows. 

• One extreme case shows an outflow with an estimated 
mass-loading factor of 7 and a momentum injection rate 
of 34 L/c. This may indicate that the observed outflow 
in this case is an overestimate, which is not unlikely 
given the large uncertainties in the observed properties. 
Otherwise, this may hint to a stronger driving mecha- 
nism. One possible way to obtain higher efficiencies is if 
the supernovae explode in extremely low density envi- 
ronments generated by the other feedback mechanism. 
Another possibility is that this clump is just now ending 
its star formation, and therefore the present measured 
SFR is smaller than the value that prevailed at the time 
most of the outflow was launched. 

We conclude that stellar feedback is expected to 
produce steady massive outflows from the high-z gi- 
ant clumps, with mass- loading factors of order unity 
and momentum injection rate efficiencies of a few. This 
is consistent with 5 of the 6-7 observed giant clumps 
where outflows were observed so far, with one or two 
exceptions in which the estimated outflows are appar- 
ently stronger. Despite the intense outflows, which in- 
dicate gas depletion times of several hundred Myr, the 
clumps arc not expected to disrupt by this process. In- 
stead, they are expected to migrate to the disc cen- 
tre on a somewhat shorter ti mescale, roughly a disc 
orbit a l time or about 250 Mvr ([Dekel. Sari fc Ceverind 
I2009t ICeverino. Dekel fc BournaudI |2010[) . and during 
this process they are expected to more than double 
in mass by accretion from the disc. One therefore ex- 
pects the population of in-situ clumps to show system- 
atic variations in their properties as a function of ra- 
dius in the disc, in the form of declining mass, stel- 
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lar age and metallicity, and in creasing gas fraction 
and specific star-forn iation rate (jCeverino et al.ll20f2l: 
iMandelker et al.ll2013( ) . One also expects that the clump 
migration, combined with the VDI-driven inter-clump 
gas in the discs, is an efficient mechanism for form- 
ing compact spheroids fPekel. Sar fc Ceverind [20091 : 

His 



ICeverino. Dekel fc BournaudI 120101: iDekel et al.M2013[ ). 
and pr oviding fuel for black- hole growth and AGN ac- 
tivity (jBournaud et al.ll201lD . 
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